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A set of quantum hydrodynamic equations are derived from the moments of the electrostatic mean-field Wigner kinetic equation. 
No assumptions are made on the particular local equilibrium or on the statistical ensemble wave functions. Quantum diffraction 
effects appear explicitly only in the transport equation for the heat flux triad, which is the third-order moment of the Wigner pseudo- 
distribution. The general linear dispersion relation is derived, from which a quantum modified Bohm-Gross relation is recovered 
in the long wave-length limit. Nonlinear, traveling wave solutions are numerically found in the one-dimensional case. The results 
shed light on the relation between quantum kinetic theory, the Bohm-de Broglie-Madelung eikonal approach, and quantum fluid 
transport around given equilibrium distribution functions. 
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1. Introduction 



Quantum effects are relevant for macroscopic charged par- 
ticle systems when the de Broglie wave-length is comparable 
^ to the mean inter-particle distance or to the size of the sys- 
^ tem, when the thermodynamic temperature is comparable to the 
Fermi temperature or for strong magnetic fields, when spin ef- 
fects are crucial. In diverse fields such as plasmas [1], semicon- 
i—{ iductors J2l or dissipative quantum models [3], quantum hydro- 
>- dynamic equations are important tools due to the relative sim- 
\Q plicity in comparison to kinetic theories and the direct physical 
,— ' interpretation of the various quantities involved. Popular ways 
to include quantum effects in the fluid equations are through 
* the so-called Bohm potential 0, Eh or by means of a quantum 
spin force |6J], among other generalized forms |7y. In these ap- 
proaches, some working hypothesis on the underlying quantum 
statistical ensemble wave functions or the local equilibrium is 
made. However, there are still issues regarding the relation be- 
tween the above different approaches that has not been clarified 
in the literature. 

In the present work, we discuss the moment hierarchy equa- 
tions derived from the electrostatic Wigner equation, which is 
the quantum counterpart of the Vlasov equation. The moments 
approach is traditional in classical kinetic theory J8(]. Without 
any further constraints or assumptions, a quantum term is ex- 
plicitly found only at the transport equation for the third-order 
moment of the Wigner function, with no simplifying assump- 
tions except that we close the system disregarding higher-order 
moments. The use of moment hierarchy quantum fluid equa- 
tions is well-known in semiconductor community [91 LlO(] but, 
to the best of our knowledge, was not applied before to the de- 
scription of electrostatic waves in quantum plasma. The reason 
for this is that in semiconductor devices there is the presence of 
a doping profile (an inhomogeneous ionic background) as well 
as an external heterojunction potential, which makes the analy- 
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sis intrinsically nonlinear ab initio. Here we assume a fixed ho- 
mogeneous ionic background and a mean-field electrostatic po- 
tential, so as to obtain a generalized linear dispersion relation, 
Eq. ( TTBT l below. A quantum modified Bohm-Gross dispersion 
relation is recovered in the long wave-length limit. The non- 
linear regimes are numerically analyzed in the one-dimensional 
case, where traveling-wave solutions are accessible. 

Furthermore, our study is also motivated as a remark on 
some recent publications II ill , in which the path from kinetic 
to fluid quantum models have been reversed. They start from 
the quantum fluid equations with a Bohm potential term, as 
derived from the Madelung decomposition of the one-particle 
wave function, and then insert this back as the momentum change 
into the Vlasov equation. However, these equations of motion 
contain the distribution function itself through the averaging 
procedure, and it is doubtful if they can be used as a basis for 
kinetic theory. It can be verified that the resulting dispersion 
relation disagrees with the Bohm-Pines dispersion relation Il2ll 
from quantum kinetic theory, except for zero-temperature equi- 
libria. Since kinetic theories are more accurate than fluid the- 
ories, these Bohmian-force methodologies are therefore con- 
troversial. See also Ref. 11 3H for related criticism. The sub- 
ject bares close connection to the present work on the moment 
quantum hierarchy, but will be more thoroughly discussed in a 
separate work. 

2. Governing equations 

Given the one-particle Schrodinger equation 



ih 



dt 2 m 



V 2 i// + e4>i// = 0, 



(1) 
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with the electrostatic potential <p, we can define the Wigner 
quasi-distribution function according to 

/(v, r) = (^j j ds exp[/mv • s/»] f(r + s/2)#- - s/2). (2) 

Then, from the Schrodinger equation, we obtain the quantum 
kinetic equation 



df 
dt 



+ V-V/ + 



j d\' K(V 



■v,r)/(v',r) = 0, 



(3) 



for the Wigner function, coupled to Poisson's equation for the 
mean-field potential, 



so 



If 



d\ f(\, r) - no 



(4) 



where no is a fixed homogeneous ionic background and K(\' - 
v, r) is defined by 

i e I tn \^ C 
K(\' - v, r) = — — — lis exp [i m (v - v') ■ s/h] 
n \2nn! J 

x[0(r + s/2)-0(r-s/2)] . (5) 

For brevity, the time-dependence of the various quantities is 
omitted. The Wigner equation (0 holds equally well in the 
case of mixed states. 

It is interesting to note that a free wave function, given by 
a Gaussian wave packet, experience a dispersive spreading (as 
expected), while the corresponding Wigner function does not 
spread for a fixed value of v 01411 . This can also be seen di- 
rectly from the structure of Eqs. (Q3 and (0. For instance, in 
the one-dimensional free-particle case and for a Gaussian ini- 
tial state if/(t = 0) = (^Jno-y 112 exp[-jt 2 /(2<x 2 )], in terms of 
a variance <x, the Wigner function can be expressed via / 
exp [-(x - v t) 2 - v 2 ]. Here the rescaled variables / = (nH/m) /, 
x = x/o~, v = mvcr/Ti, and t - ht/(mcr 2 ) are employed. The 
corresponding contour-plot graphics are shown in Fig. 1. 

2.1. Moment equations 

To obtain macroscopic equations, let us introduce the mo- 
ments 

n = J dvf, (6) 
«u = Jd\ f v , (7) 
Pij = m ^J"dvfvi vj - n UjU^j , (8) 
Qijk = m Jd\ (vi-Ui)(vj-Uj)(v k -u k )f, (9) 
Rijki = m Jdv(Vi - iii)(vj - uj)(v k - u k )(vi - m;)/,(10) 

from which, in particular, we can derive a scalar pressure p = 
(1/3) Pa and a heat flux vector = (1/2) Qjp Here and in the 
following, the summation convention is applied. The Wigner 
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Figure 1: Time-evolution of the free-particle rescaled Wigner function / = 
(nti/m)f, in terms of non-dimensional variables x = x/cr,v = mvtr/Pi and 
i = Ht/(m(T 2 ). Initial wave function: i//(t = 0) = (-0rcr) -1 ' 2 exp[-,r 2 /(2<T 2 )]. 
Upper, left: t = 0; upper, right: t = 2; bottom, left: t = 4; bottom, right: 1=6. 



equation then gives the following macroscopic equations, 



Dn 

= -nv ■ u , 

Dt 

Dm djP u e 
— = + - di4> , 

D t mn m 



(11) 
(12) 



DPij 



Dt 

D Qi jk 1 



P Hi d k uj ) -P ij V-u-d k Q i j l 



(13) 



,, p (U d ' p k)i ~ Qi(ijd'uk) ~ QijkV ■ u 

Dt mn 

eh 2 n 

ijkl 



Am 2 ijk 



(14) 



where d,- = 6ijd^ = d/d r,- and the material derivative is D/Dt = 
d/dt + u-V. The calculation assumes e.g. decaying or periodic 
boundary conditions and uses the symmetry properties of Pjj, 
Qij k and /?y# under permutation of indices. Finally, in Eqs. 
([T3l- (fT4l i. the round brackets denote symmetrization, where we 
use a minimal sum over permutations of free indicices needed to 
get symmetric tensors. Thus, for example, with this convention 
P k (id k Uj) is defined as P k ^ d k Uj) = P k i d k Uj+Pj k d j Ui+Pijd'u k . 

We here remark that a) there is no assumption on the par- 
ticular equilibrium Wigner function. This is a difference in 
comparison with some previous approaches [4] relying on the 
first-order quantum correction to Maxwell-Boltzmann equilib- 
ria Therefore, the model is not semi-classical and is not 
restricted to classical statistics; b) the explicit dependence on 
Planck's constant appear only when the heat flux triad trans- 
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port equation is considered. In addition, the quantum contribu- 
tion disappears in Eq. ( TBI when the scalar potential is absent. 
This is similar to Gardner's approach and in contrast to the 
method of Ref. [5], where quantum effects modeled by a Bohm 
potential appear already at the momentum transport equation, 
via the pressure dyad Py and the associated Madelung decom- 
position of the quantum statistical ensemble wave functions. 
Here the usual quantum force is replaced by the third-order 
derivative of the scalar potential term in Eq. (TT~4-b . 

It is to be expected from the very beginning that Planck's 
constant would not appear through the moments of the Wigner 
equation when there is no electric field, because in this case the 
Wigner equation reduces to the free-particle Vlasov equation 
(see the discussion above and Fig. 1) and the initial conditions 
on the Wigner function determines the quantum aspects of its 
evolution. To conclude, in the present context there are two 
sources of quantum terms, one via the explicit dependence of 
the Wigner equation in n and the other via the definition of the 
Wigner function associated with the proper quantum statistical 
ensemble. 

2.2. Closure and dispersion relation 

In this work we choose the simplest way to achieve clo- 
sure of the system (fTTT)-([T4l). neglecting the contribution from 
the fourth-order moment Ryu- We take into account Poisson's 
equation and linearize around the homogeneous equilibrium n = 
no, u = 0, 



[Pij] = n K B T Q ±(x ® x + y ® y) + r ||Z ® z 



(15) 



Qijk — 0, and = 0, where the equilibrium temperatures per- 
pendicular and parallel to the wave propagation Tq ± and To\\ can 
be unequal. Here, Kb is Boltzmann's constant and plane wave 
perturbations proportional to exp(ikz - icot) are assumed, with- 
out loss of generality. It follows that 



2 W P 
U> = 



1 + 



1 



12 k b T »k 2 h 2 k 



2lA \ 



m 2 oj 2 t 



1/2 



(16) 



where u> p = (e no/meo) 1 is the plasma frequency and the 
equilibrium temperature Tqj. perpendicular to the wave prop- 
agation does not contribute. In the particular case of small 
wave-vector and quantum effects, Eq. ( [Tol l reduces to the usual 
quantum Langmuir dispersion relation a> 2 = a) 2 +3 Kb Tq k 2 /m+ 
ti 2 k 4 /4 m 2 . However, even in the formal classical limit (H = 0), 
Eq. ([TBI yields the Bohm-Gross dispersion relation only for 
long wave-lengths. Finally, despite the appearance, the new dis- 
persion relation is not restricted to Maxwell-Boltzmann equilib- 
ria. For instance, in the case of a zero-temperature, degenerate 
electron gas, the only basic change would be the substitution of 
the parameter Toy by the Fermi temperature. 

Contrarily to the habitual usage, a scalar pressure dyad Py = 
p 6ij is not a valid assumption even in our electrostatic case. In- 
deed, linearization of Eqs. (fT2l-(fT4li with Py = Poy + £SP,j, 
<p = e6<p, £ <K 1, gives 



SPy = - 



e&bk 2 ( n n 2 k 2 6 iz 6 jz \ 
-^y [Poy + P 0(iz 6 jz) + I (17) 



m oj- 



4 m 



as the first-order perturbation of the pressure dyad, assuming 
kj = k 5j Z without loss of generality. Clearly the wave propaga- 
tion itself is a source of anisotropy, even for isotropic equilib- 
ria and/or purely classical plasma. The result dTTb also follows 
from the kinetic theory in the long wave-length limit. Further- 
more, as apparent from Eq. (fT7b , it is legitimate to postulate 



[Py] = hk b T ± (t ® x + f ® f) + T\\Z ® z 



(18) 



where the perpendicular and parallel temperatures T ± and T\\ in 
general are different. 

We note that the Bohm-Gross dispersion relation is recov- 
ered from Eqs. (fTT1 i-(fT3b in the adiabatic and classic case. In- 
deed, setting Qyit = and linearizing, the result is to 2 = co 2 + 
3 (k b 7on/m) k 2 , where only the component = no k b Tq\\ of 
the equilibrium pressure dyad contributes. On the other hand, 
insisting on an isotropic pressure dyad Py = pdij and tak- 
ing Qijk = 0, and combining the continuity equation with Eq. 
( fT3l we get p = noK B To(n/no) 7 , where y = 5/3, implying 
the dispersion relation a> 2 = a> 2 + y(KBTo/m)k 2 . To prop- 
erly recover the Bohm-Gross dispersion relation in an adiabatic, 
scalar pressure fluid theory retaining only up to the first or- 
der moment, there is the need of a phenomenological adiabatic 
exponent y = 3, reflecting the fact that plane wave propaga- 
tion is essentially a one-dimensional phenomena 11611 . Alter- 
native moment hierarchy formulations II 1 Oil , closed at the tem- 
perature (basically the trace of the second-order moment of the 
Wigner function) evolution equation, can be shown to result in 
cj 2 = c^ + (5/3) (K B To/m) k 2 +ti 2 k 4 /(12 m 2 ), which goes neither 
to the Bohm-Gross nor Bohm-Pines dispersion relations in the 
classical or zero-temperature limits, respectively. In contrast, as 
shown here, the quantum modified Bohm-Gross dispersion re- 
lation is a natural consequence from third-order moment theory, 
in the long wave-length limit of Eq. ( fTol l. 

2.3. Dynamics 

As an example of the dynamics of the third order hierarchy 
model, we look at the strictly one-dimensional case, such that 
(di = 8 ix d x , Ui = u6 ix , Py = p5 ix 6 jx , Qy k = QS ix S jx S kx ), 



and 



h + n d x u — , 

d x p e d x <p 

u = 1 , 

mn m 

p = -3pd x u -d x Q, 

. 3 p d x p eh 2 n d|g> 

Q = -7—2 4Qd x it, 

m n 4 m z 



d x <f> = — (n - n ) , 
so 



(19) 
(20) 
(21) 

(22) 
(23) 



where the dot represents the convective derivative and disre- 
garding the fourth-order moment. The dispersion relation in 
Eq. <[16J follows from Eqs. (HMS). 

Assuming traveling solutions with all quantities depending 
only on the variable £ = x— vt, where v is a fixed parameter, the 



following nonlinear system of ordinary differential equations is 
obtained, 



{u — v) u 
(u - v) p' 

(H-V)G' 



mn m 
-3pu' - Q' , 
3 p p' eh 2 n< 



4 m 2 



-4Qu', 



so 



-(n - n ) . 



(24) 
(25) 

(26) 

(27) 



where the prime denotes derivative with respect to Finally, 
the continuity equation can be integrated to hqUq = n(u - v) 
etc.. in terms of a reference velocity uq which we assume to be 
nonzero to exclude trivial cases. 

Eliminating n through the continuity equation, it can be 
shown that the system of Eqs. (I24l >-(|27]> admit linearly stable 
oscillations around the equilibrium u = uq + v, p = po,Q = 
0, <p — <p' = 0, provided the inequality hu)pl(mufy < 2 is sat- 
isfied. Figures 2 and 3 show typical oscillations in this case. 
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Figure 2: Particle density and velocity field oscillations from the system of Eqs. 



24H27t . where f = x - vt, Ko + and for ti w p /(m iify = 1. Initial conditions 



such that «(0) ; 
O,0'(O) = 0. 
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Figure 3: Pressure, heat flux and electric field oscillations from the system 
of Eqs. )24H27t . where f = x — vt, uq + 0, and for ha) p l(mu^) = 1. Initial 
conditions such that «(0) = (2/3) no, u(0)- v = (3/2) »o,p(0) = m«o«Q,2(U) = 
0,^(0) = 0,f (0) = 0. 



of a Madelung decomposition of the underlying quantum statis- 
tical ensemble wave functions or assumptions on the local equi- 
librium configuration. A generalized dispersion relation for lin- 
ear waves is derived, from which the quantum modified Bohm- 
Gross dispersion relation is recovered in the long wave-length 
limit. For closure, fourth-order moments were discarded. More 
sophisticated closure schemes |l7ll designed to reproduce some 
of the results from kinetic theory, as well as the inclusion of 
spin effects, are postponed to future considerations. 



3. Conclusions 

To conclude, a higher order moment quantum hydrodynamic 
model for quantum plasmas was derived, starting from the elec- 
trostatic Wigner equation. Quantum effects appear explicitly 
going up to the third-order moments hierarchy, without the need 
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